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Abstract. In future measurements of the dilepton (Z/7*) transverse momentum, Qt, at both the Tevatron 
and LHC, the achievable bin widths and the ultimate precision of the measurements will be limited by 
experimental resolution rather than by the available event statistics. In a recent paper the variable ar, 
which corresponds to the component of Qt that is transverse to the dilepton thrust axis, has been studied 
in this regard. In the region, Qt < 30 GeV, ar has been shown to be less susceptible to experimental 
resolution and efficiency effects than the Qt- Extending over all Qt, we now demonstrate that dividing 
(it (or Qt) by the measured dilepton invariant mass further improves the resolution. In addition, we 
propose a new variable, </>^, that is determined exclusively from the measured lepton directions; this is 
even more precisely determined experimentally than the above variables and is similarly sensitive to the 
Qt- The greater precision achievable using such variables will enable more stringent tests of QCD and 
tighter constraints on Monte Carlo event generator tunes. 

PACS. 12.38.Qk - 13.85. Qk - 14.70.Hp 



1 Introduction 

The production of Z/j* — > e + e~ and Z/j* — > fi + fi~ at 
hadron colliders provides an ideal testing ground for the 
predictions of QCD, due to the colourless and relatively 
background free final state. The dilepton transverse mo- 
mentum, Qt, distribution probes QCD radiation in the 
initial state. At high values of Qt (Qt > 30 GeV, say), 
the fixed order perturbative calculations now available at 
NNLO [3(2] are expected to yield accurate predictions. 
At low Qt, soft gluon resummation techniques are re- 
quired [3J, with additional non-perturbative form factors 
determined in global fits to experimental data such as 
in [4]. Various event generators are also available [3[6j[7l 
[5] , matching tree level matrix elements to parton showers 
tuned to data. Validation and tuning (form factors and 
parton showers) of these models require comparison with 
experimental data that have been corrected for detector 
resolution and efficiency effects. Improved understanding 
of these production models will increase sensitivity to new 
physics signals at hadron colliders. 

The Qt distribution has been measured at the Fermi- 
lab Tevatron, by the CDF and D0 [TO1[TT1[T2] Collabo- 
rations. The most recent of the above measurements [TTJ 
IT2] used approximately 1 fb" 1 of data. Although this rep- 
resents only about one tenth of the anticipated final Teva- 
tron data set, the precision of these measurements was 
already limited by experimental systematic uncertainties 



in the corrections for event selection efficiencies and un- 
folding of lepton momentum resolution. In order to un- 
fold measured distributions for experimental resolution it 
is important that the chosen bin widths are not too small 
compared to the experimental resolution. In the low Qt 
region in [111112] , the minimum bin sizes were determined 
by experimental resolution rather than the available data 
statistics. The final Tevatron data set will be an order of 
magnitude larger than that analysed in [Tl"l[T2"] . New ideas 
are therefore needed in order to exploit fully the data for 
studying the physics of boson Qt- 

The olt variable was introduced in |13j and was pro- 
posed as a novel variable for studying the Qt in 14 . Fig- 
ure [T] illustrates this and other relevant variables defined 
below. Events with A(j> > tt/2, where A(f> is the azimuthal 
opening angle of the lepton pair, correspond to approx- 
imately 99% of the total cross section. For such events 
the Qt is split into two components with respect to an 



event axis defined as, i = [pt 



(i) 



Pt 



(2) 



)/\PT 



(1). 



Pt 



(2) 



where Pt^ and Pt^ are the lepton momentum vectors 
in the plane transverse to the beam direction. The com- 
ponent transverse to the event axis is denoted by ay and 
the aligned component is denoted by a^. For events with 
A(j) < tt/2 this decomposition is not useful and or and 
cil are defined as being equal to Qt- clt is less susceptible 
than Qt to the lepton pt resolution. In addition, the ef- 
ficiencies of selection cuts on lepton isolation and pt are 
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shown in [2] to be less correlated with or than Qt- For 
studying the low Qt (non-perturbative) region, or is thus 
a more powerful variable than Qt- The cit distribution has 
subsequently been calculated to NLL accuracy using soft 
gluon resummation techniques |15j . 




Fig. 1. Graphical illustration in the plane transverse to the 
beam direction of the variables defined in the text and used to 
analyse dilepton transverse momentum distributions in hadron 
colliders. 

A recent paper [16] has discussed the idea of using A(f>, 
as an analysing variable that is sensitive to the physics 
of Qt, and insusceptible to lepton momentum uncertain- 
tieqj. Whilst A<\> is primarily sensitive to the same compo- 
nent of Qt as ar, the translation from q,t to A(f> depends 
on the scattering angle 9* of the leptons relative to the 
beam direction in the dilepton rest frame. Thus, Acj> is 
less directly related to Qt than <zt- As a result, A<j> has 
somewhat smaller statistical sensitivity to the underlying 
physics than aT- 

In this paper, we discuss two further ideas to improve 
experimental precision, whilst maintaining (Qt) physics 
sensitivity: 

— Dividing ax, and Qt by the dilepton invariant mass, 
Q, thus further reducing the effects of lepton pt res- 
olution, and almost totally cancelling lepton pt scale 
calibration uncertainties. 

— Correcting Acf> on an event-by-event basis for the scat- 
tering angle, 9* , thus improving the sensitivity to Qt- 

An overview of the rest of this paper is as follows. In 
Sect.[5]we give an approximate analytic motivation for the 
idea of dividing ot (and Qt) by Q in order to produce 
variables with very substantially improved experimental 
resolution. In Sect.[3]we discuss the idea of correcting A<j) 
on an event-by-event basis for the scattering angle, 9*, 
thus improving the sensitivity to Qt- In addition, we pro- 
pose a new variable, cos(#*), which provides a measure 
of the scattering angle that is based entirely on the mea- 
sured track directions and is thus extremely well measured 
experimentally. In Sect. [4] we describe the simple param- 
eterised detector simulation we employ in our MC stud- 

1 We note that the expected distribution of A(f> does have a 
small residual sensitivity to the lepton pr measurement. This 
arises from the cut on Q in the event sample selection, which 
is affected by the lepton pr scale and resolution. 



ies. In Sect. [SHHlwe present the results of our MC studies 
of the performance of the various candidate variables in 
terms of their experimental resolution and their sensitiv- 
ity to the underlying physics. In Sect. [TU] we present some 
conclusions, including our recommendations for the best 
variables to use in experimental studies of the transverse 
momentum of dilepton pairs produced at hadron colliders. 



2 Mass ratios of a T and Q T 

For A<f> w 7T, aT is given by the approximate formula 



Pt +Pt 

and thus the fractional change in ay with respect to a 
variation in, say, p/p is given by 

Aar _ P ( t Ap ( t 
The dilepton invariant mass is given by 

Q = yVV^C 1 - co s 

where p^ and p^ are the lepton momenta and A9 is 
the angle between the two leptons. Thus, the fractional 
change in mass with respect to a variation in p^ is given 

by 

AQ _ 1 Ap^ 
~~Q ~ 2 p« ' 

Since track angles are extremely well measured it can be 
taken to a very good approximation that 

~pY~~p^' 

The fractional change in aT /Q with respect to a variation 
in p^ is thus given by 

A(a T /Q) _ Aar _ AQ _ ( l\ Apfi 

(clt/Q) ~ a T Q ~ [pW+pW ~ 2) 

Thus the variations with pj? in or and Q partially 
cancel in the ratio, rendering ay/Q less susceptible to the 
effects of lepton pt resolution than ay. In particular, in the 

region of low Qt then p^p « p^ and thus A(o.t/Q) ~ 0. 
Similarly, the quantity Qt/Q is less susceptible to the 
effects of lepton pt resolution than Qt- 

A simple example of an uncertainty in the lepton pt 
scale calibration is to consider the pt of all leptons to be 
multiplied by a constant factor. It can be seen trivially 
that in this case aT, Qt and Q are all multiplied by the 
same factor and that the measured clt/Q and Qt/Q are 
unaffected by such a scale uncertainty. 
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3 Correcting A<p for the scattering angle 

The azimuthal opening angle between the two leptons, 
A(f>, is primarily sensitive to the same component of Qt as 
cit, and is based only on the well measured lepton angles. 
However, at fixed <it/Q, A<f> depends on the scattering 
angle 9* of the leptons relative to the beam direction in 
the dilepton rest frame. For convenience, we define the 
acoplanarity angle, acO p, as acop = n - A(f>. For pf> w 

(2) 

p T it can be fairly easily shown that 

a T /Q~ tan(^ acop /2) sin(0*). 

This suggests that the variable 

cf)* = tan(^ acop /2)sin(0*) 

may be an appropriate alternative quantity with which to 
study Qt- 

In the analysis of hadron-hadron collisions, 9* is com- 
monly evaluated in the Collins Soper frame [17] . However, 
9q S requires knowledge of the lepton momenta and is thus 
susceptible to the effects of lepton momentum resolution. 
Motivated by the desire to obtain a measure of the scat- 
tering angle that is based entirely on the measured track 
directions (since this will give the best experimental res- 
olution) we propose here an alternative definition of 9*. 
We apply a Lorentz boost along the beam direction such 
that the two leptons are back-to-back in the r-9 plane. 

This Lorentz boost corresponds to ft = tanh ij 3 — ) and 
yields the resuli0 



the particle level momenta: 5(1 /px) — 3 x 10~ 3 (1/GeV) 
for muons, which are measured in the tracking detectors; 
Sp/p — 0.4(p/po) -1 / 2 with po = 1 GeV for electrons, 
which are measured in the calorimeter. In addition, the 
particle angles are smeared, assuming Gaussian resolu- 
tions of 0.3 x 1CT 3 rad for (j> and 1.4 x 10~ 3 for r\. These 
energy, momentum, and angular resolutions roughly cor- 
respond to those in the D0 detector [2Tj . 

Events are accepted for further analysis if: 70 < Q < 
110 GeV and both leptons satisfy the requirements pt > 
15 GeV and \rj\ < 2. These selection cuts are made at 
particle level, unless otherwise stated. 



5 Scaling factors 

In the following sections, we compare the experimental 
resolution and physics sensitivity of the various candidate 
variables. In particular, we compare the variation of the 
resolution for each variable as a function of that variable. 
This comparison is facilitated by ensuring that all dis- 
tributions have approximately the same scale and shape. 
Compared to Qt or Qt/Q, all other variables are on av- 
erage a factor y/2 smaller (since ay and measure one 
component of Qt)- A simple multiplication by Mz (= 
91.19 GeV [32]) corrects for the average Q _1 factor in the 
mass ratio and angular variables and conveniently ensures 
that all variables have the same units (GeV). Finally, the 
mean value of sin(#*) is around ~ 0.85, and tan(</> acop /2) 
is scaled by this additional factor. The above factors are 
summarised in Table [TJ 



tanh 



7} 



where rf~ and i] + are the pseudorapidities of the negatively 
and positively charge lepton, respectively. 
We consider two candidate variables 



^cs 



tanf 



3 /2)sin(0* s ) 



4>* v = tan(0 acop /2) sin(0*) 

for further evaluation in terms of their experimental reso- 
lution and physics sensitivity. 



variable scaling factor 



Qt 
Qt/Q 

ax 
ar/Q 

ah/Q 
tan(0 acop /2) 

<^cs 



1 

M z 

V2 
V2M Z 

V2 
V2M Z 
0.85y/2M z 
V2M Z 
V2M Z 



Table 1. Scaling factors for different candidate variables. 



4 Simple parameterised detector simulation 

Monte Carlo events are generated using pythia 18 , for 
the process pp — > Z/j*, in the e + e~ and decay 
channels, and re-weighted in dilepton Qt and rapidity, 
y, to match the higher order predictions of resbos [TH] 
as in [TJ]. Electrons and muons are defined at "particle 
level" according to the prescription in [20], and at "de- 
tector level" by applying simple Gaussian smearing to 

2 The lepton pseudorapidity, r), is defined as r) — 
— ln[tan(!)], where 6 is the polar angle with respect to the 
beam direction, in the laboratory frame. 



6 Experimental resolution for dilepton 
scattering angle 

Figure[5]shows the experimental resolution of cos (9 £ s ) and 
cos(0*) in our simulation of dimuon events. The upper row 
of Figure [5] shows events that satisfy 70 < Q < 110 GeV; it 
demonstrates that cos(6>*) is significantly better measured 
experimentally than cos(0J g ). This is because cos(9*) is 
evaluated using only angular measurements, which are 



4 



A. Banfi et al.: Optimisation of variables for studying dilepton transverse momentum distributions 



more precise than the momentum measurements included 



in the determination of cos( 



'cs 



is used in the definition of <j>* = 



The variable cos( 
tan((/) acop /2) sin(6>*) in Section [3] above. As an aside, we 
note in addition that a precise determination of the dilep- 
ton centre of mass scattering angle that is free from the 
effects of lepton momentum resolution can also find ap- 
plication in the determination of the forward-backward 
charge asymmetry of dilepton production at hadron col- 
liders. The experimental resolution in cos(#£ g ) becomes 
particularly significant in the dimuon channel for very high 
values of Q for which the lepton momentum resolution is 
poorest. This is illustrated in the lower row of Figure [21 
which shows the experimental resolution of cos(6*£. s ) and 
cos(6»*) in events that satisfy 500 < Q < 600 GeV. The 
advantage of using cos(#*) for high mass events is even 
larger than was the case for 70 < Q < 110 GeV. 



7 Experimental resolution for variables 
related to dilepton Q T 

Figure [3] compares separately for our dimuon and dielec- 
tron simulations, the mean resolution, \Sx\, of various can- 
didate variables, function of (particle level) x. All 
variables are scaled by the factors in Table [TJ 
The following observations are made: 

— o-t/Q is significantly better measured than ax, over 
the entire range. 

— Similarly, Qt /Q is significantly better measured than 
Qr- 

— Over the full range, ot and ax/Q perform better than 
Qt and Qt/Q respectively. 

— Compared to cit/Q, fics nas significantly better reso- 
lution, and 0* better still. 

— The most precisely measured variable is tan(<^ acop /2), 
since it is determined only from the azimuthal angles 
of the leptons, whereas the uncertainty on </>* includes 
also the uncertainties on the measured lepton pseudo- 
rapidities. 

Since the discussion in Section [5] is only approximate, 
we have investigated empirically various other possible 
scalings of cit with Q (with the appropriate scale factor 
applied, as above) . These are illustrated in Figure 2J As 
expected, it can be seen that when compared to the other 
variables considered in Figure El ar/Q has the best exper- 
imental resolution for all Qt and irrespective of whether 
(a) tracker-like or (b) calorimeter-like resolution in the 
lepton momenta is simulated. 



8 Event selection efficiency 

As discussed in [TJ], the efficiencies of selection cuts on 
lepton isolation and pt (for ot < P™') are l ess correlated 
with a,T than Qt- We have verified that the correlations 
of these selection efficiencies with cit/Q and Qt/Q are 
essentially identical to those with q,t and Qt respectively. 
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(a) 70 < Q < 110 GeV 
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(b) 70 < Q < 110 GeV 
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5 cos(e*) 

(d) 500 < Q < 600 GeV 

Fig. 2. Distributions of experimental resolutions the dilepton 
centre of mass scattering angle for events satisfying 70 < Q < 
110 GeV (upper row) and 500 < Q < 600 GeV (lower row). 
Figures (a) and (c) show cos(#£s)- Figures (b) and (d) show 
cos(6<*). 
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Particle level x (GeV) 

(a) Tracker 
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10" 2 50 100 

Particle level x (GeV) 

(b) Calorimeter 

x = 

Q T Q T /Q (x M z ) 

a T (xV2) a T /Q (xV^IVy 

a L (x\/2) a L /Q (xV^IVy 

^(xN/ay ^(xVsm,) 



tan((|) ac 12) (x 0.85 xV^A/y 



Fig. 3. The mean resolution of each candidate variable, as a 
function of that variable (scaled by the factors described in the 
text). Results are presented both for (a) tracker- like and (b) 
calorimeter-like resolution in the lepton momenta. 
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X = 




a T C^) 


-- a T /Q (xV^ivy 


- a T AfQ (xV^) — 


— a T /Q 2 (x\l2M z ) 


a T /ln(Q/Q ) (xV2ln(M IQ Q )) 





Fig. 4. The mean resolution of the variables ar, ar/Q, 
a T /Q 1/2 , clt/HQ/Qo), and a T /Q 2 , as a function of that vari- 
able (scaled by the factors described in the text). Results are 
presented both for (a) tracker-like and (b) calorimeter-like res- 
olution in the lepton momenta. 



6 



A. Banfi et al.: Optimisation of variables for studying dilepton transverse momentum distributions 



A(j> is primarily sensitive to the ax component of the Qt- 
The variables, </>*, and 4>csi are verified to exhibit the 
same benefits as ax compared to Qt, in terms of efficiency 
dependence. 



9 Sensitivity to the physics 

Figure O shows the particle level, normalised distributions 
of a T , M z a T /Q and ln(M z /Q )a T / HQ/Qo) (with Q = 
1 GeV) for three ranges of Q. We see that ar has a mild 
dependence on Q, while dividing by Q over corrects this 
dependence. In this respect, we observe that the distribu- 
tion in ax / ]n(Q / Qa) has a smallest dependence on Q, as 
might be expected from [3J. 

Experimental measurements of Z boson production 
are typically made over a fairly wide bin in Q (e.g., 70- 
110 GeV). One potential concern with measurements of 
clt/Q and Qt /Q is that the increased correlation with Q 
demonstrated in Figure [5] might degrade the sensitivity to 
the underlying physics. Since (j)* behaves approximately 
as (/)* rs cit/Q, a similar degradation in the physics sen- 
sitivity of (/>* may similarly be expected. In this respect, 
Or/ ln(Q/Qo) is a more suitable variable than clt/Q for 
studying the boson Qt distribution over a wide range in 
Q. However it has poorer experimental resolution, as was 
demonstrated in Figure |U 

Using a similar procedure to that described in [T4l . 
we study the sensitivity to the underlying physics of the 
different candidate variables. We run pseudo-experiments 
to fit for the value of the parameter <?2, which determines 
the width of the low Qt region in resbos. Events must 
meet the requirements described in Section [4j except that 
these are applied at detector level rather than particle 
level. Of these events, 1M are assigned as pseudo-data and 
the remaining events are used to build g2 templates. All 
variables are scaled by the factors listed in Table [TJ such 
that the same binning (30 equal width bins in the range 
0-30 GeV) can be used. A minimum \ 2 fit determines the 
statistical sensitivity of each variable to the value of <?2 • 

Whilst in [TT] we were mostly interested in the low Qt 
region, the ideas proposed in this paper also improve ex- 
perimental resolution at higher Qt (see Section [7]) ■ Thus 
it is of interest to study the physics sensitivity also in 
this region. We similarly fit for a parameter Kq t , which 
weights events with (particle level) Qt > 25 GeV, by 
Kq t {Qt — 25), and approximately represents the differ- 
ences between predictions at NLO and NNLO discussed 
in [TT] . Again, after applying the appropriate scaling fac- 
tors from Table [TJ the same binning can be used for each 
variable^. 

The results of the fits to 32 and Kq t are presented in 
Tables [U and respectively. Results are given separately 
for particle-level (dimuon) and detector-level (tracker and 
calorimeter). For both o,t/Q and Qt/Q, the statistical 

3 The first bin is of width 5 GeV (with lower edge at zero) 
and each consecutive bin is 2 GeV wider than the last. Ten 
such bins give an upper edge of the last bin at 140 GeV and 
the fit includes the overflow bin from 140 GeV to 00. 
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2 4 6 8 10 12 14 16 18 20 

M z a T /Q (GeV) 




2 4 6 8 10 12 14 16 18 20 



ln(M z /Q )a T /ln(Q/Q Q ) (GeV) 

Fig. 5. Comparison of particle level distributions of ar, 
Mzar/Q and \n(Mz /Qo)ar / hi(Q/Qo) for three ranges of Q. 



sensitivities are essentially the same as those for or and 
Qt respectively. Thus the effect of the additional Q de- 
pendence is shown to be negligible. 

The approximately 5% poorer sensitivity of tan(</> acop /2), 
compared to ay/Q, demonstrates the sin(#*) ambiguity 
of the former. The additional factor sin(#*), present in <f>* 
and i^q S actually recovers the sensitivity to the same level 
as aT- In addition, the <f>* variable, which was shown in 
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variable 


particle level 


calorimeter 


tracker 


Qt 


0.65 


0.94 


1.41 


Qt/Q 


0.65 


0.94 


1.40 


(It 


1.00 


1.00 


1.00 


a T /Q 


1.00 


1.01 


1.00 


ah 


1.21 


2.35 


4.74 


tan(</> acop /2) 


1.04 


1.05 


1.04 


^cs 


1.00 


1.00 


0.99 


K 


1.00 


1.00 


0.99 



Table 2. Statistical sensitivity (in %) on the parameter (?2 from 
fits to the distributions of different of variables. For details see 
text. 



variable 


particle level 


calorimeter 


tracker 


Qt 


1.65 


1.67 


1.82 


Qt/Q 


1.66 


1.67 


1.78 


ax 


1.92 


1.92 


1.96 


a T /Q 


1.92 


1.92 


1.94 


ah 


1.98 


2.02 


2.34 


tan(</> acop /2) 


1.96 


1.96 


1.98 




1.88 


1.88 


1.90 


K 


1.87 


1.87 


1.92 



Table 3. Statistical uncertainty (in %) on the parameter Kq t 
(as defined in the text) from fits to the distributions of different 
of variables. For details see text. 



Section [7] to have the best experimental resolution (ex- 
cept for tan((/) acop /2)), performs similarly to 4>q S in terms 
of physics sensitivity. 

Of course, the results presented in Tables [5] and [3] rep- 
resent only the statistical sensitivity of the considered vari- 
ables when compared to Qt- As discussed in [T3], the sys- 
tematic uncertainties associated with modelling of detec- 
tor resolution and efficiency will be significantly smaller 
using variables that are better measured and less corre- 
lated with event selection efficiencies than is the case for 
Qt- 



10 Summary and conclusions 

Measurements of dilepton (if/7*) transverse momentum, 
Qt, distributions are crucial for improving models of vec- 
tor boson production at hadron colliders. The precision of 
future measurements of the Qt distribution at the Teva- 
tron and LHC will be totally limited by uncertainties in 
correcting for detector resolution and efficiency, and the 
minimum bin sizes will be limited by resolution. 

In |14j an alternative variable, ax, was demonstrated 
to be significantly less susceptible to such detector effects 
than Qt- We have shown in this article that the experi- 
mental resolution of ot can be further improved by tak- 
ing the ratio to the measured dilepton invariant mass, Q. 
Similarly, we have demonstrated that the variable Qt/Q 
is experimentally more precisely measured than Qt- No 
obvious disadvantages of the variable ax/Q (Qt/Q) as 
compared to (Qt) are found when studying the effi- 



ciency dependencies or the physics sensitivity (by making 
fits to parameters describing the shape of the Qt distri- 
bution) . 

The acoplanarity angle, acO p, is also sensitive to Qt- 
However, it has the approximate dependence Ot/Q ~ 
tan(</> acop /2) sin(0*), where 8* is the scattering angle of 
the leptons with respect to the beam direction in the dilep- 
ton rest frame. Thus </> ac0 p is less directly related to Qt 
than or- We show that correcting tan(</> acop /2) by factor 
sin(0*) yields a variable with essentially the same statis- 
tical sensitivity as aT/Q, whilst the experimental resolu- 
tion is significantly better. Furthermore, using an approx- 
imate rest frame determined using only angular informa- 
tion, we propose a new method to measure the scattering 
angle, 9*, with the best possible experimental resolution. 
We conclude that in the region of low Qt the variable 
(f>* = tan(0 acop /2) sm(9*) represents the optimal combi- 
nation of physics sensitivity, experimental resolution and 
immunity to experimental systematic uncertainties. 

For studying the high Qt region, the optimal variable 
is Qt/Q, which is significantly better measured than Qt, 
and has no disadvantages in terms of physics sensitivity or 
efficiency dependence. However, it is interesting to study 
correlations between the ay and ax components of Qt- 
Therefore, measurements of (f>* in the high Qt region will 
be complementary to measurements of Qt/Q- 

These further optimisations of variables used to study 
the Qt distribution will allow significantly finer binning 
and smaller unfolding corrections (and thus uncertainties) , 
enabling tighter constraints on vector boson production 
models. 
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